Conductance fluctuations and weak localization in chaotic quantum dots 
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We study the conductance statistical features of ballistic electrons flowing through a chaotic 
quantum dot. We show how the temperature affects the universal conductance fluctuations by 
analyzing the influence of dephasing and thermal smearing. This leads us to two main findings. 
First, we show that the energy correlations in the transmission, which were overlooked so far, are 
important for calculating the variance and higher moments of the conductance. Second, we show 
that there is an ambiguity in the method of determination of the dephasing rate from the size of 
the of the weak localization. We find that the dephasing times obtained at low temperatures from 
quantum dots are underestimated. 
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A very striking experimental evidence of universal sta- 
tistical behavior due to quantum coherence and complex- 
ity in electronic ballistic transport was recently reported 
by Huibers and collaborators |Q. They measured the 
conductance G through a chaotic quantum dot at small 
bias and low temperatures as a function of an applied 
magnetic field and the quantum dot shape. For such de- 
vices, where the quantum coherence length 1$ and the 
system size L are such that <^ » i, the conductance 
is expected to parametrically display mesoscopic fluctu- 
ations. 0] To characterize the latter, due to the system 
complexity, a detailed microscopic theory is neither fea- 
sible nor practical. Hence, the indicated theoretical ap- 
proach should be statistical and tailormade to give the 
experimental accessible statistical measures such as the 
conductance distribution P(G), conductance autocorre- 
lation functions, etc.. For ballistic chaotic quantum dots 
such approach is provided by the random matrix theory 
(RMT). Indeed, the agreement between the conductance 
distributions P{G) obtained in Ref. Jl|] and the corre- 
sponding stochastic theory, turned this experiment into 
a paradigm of the statistical approach. |||| 

Early experiments revealed an unexpected aspect 
to that systems, namely that even at low temperatures 
the conductance fluctuations significantly deviate from 
the predictions of the simplest random matrix models 
(@,|). More specifically, we are referring to the suppres- 
sion of the weak localization peak, which represents the 
first quantum correction to the classical picture, and to 
the conductance variance var(G). The early works were 
improved and converged to the understanding that even 
a small loss in quantum coherence J^] affects dramati- 
cally the statistical observables At the quantitative 
level, some features of the experimental data still remain 
unexplained. 

The main findings presented in this letter are two-fold. 
First, using an alternative statistical approach we explain 
the discrepancy between theory and experiment in Ref. 



jjj for var(G). This result has important consequences 
for recent predictions of var(G) in similar systems. [To| ] 
Second, we show that there is an ambiguity in the way 
the dephasing rates are extracted from the weak local- 
ization experimental data in open chaotic quantum dots 
so far. (ll]-[l3) Within our statistical model we propose 
a different method, which indicates that the dephasing 
rates quoted in the literature [jl2| are overestimated. 

The conductance G = (e 2 /h)g through a two-lead 
quantum dot is related to the transmission, and hence 
to the S matrix, by the Landauer formula 



g(E, X) = T 21 (E, X) = Y] \S ba (E, X) 



(1) 



Here g is the dimensionless conductance, T21 is the trans- 
mission of an electron scattered from the incoming lead 
1 to the outgoing lead 2, and the labels of the corre- 
sponding scattering matrix S indicate the open channels 
located at each lead. X is a generic parameter such as 
a gate voltage, which shapes the dot, or an external ap- 
plied magnetic field B. The applicability of the Landauer 
formula assumes full quantum coherent transport. 

Thermal effects modify Eq. ([!]) in different ways. First, 
and most interesting, by increasing the temperature the 
dynamics in the dot changes, making the coherent single- 
particle description of the process less realistic. The rich 
physics involved attracted a lot of attention and a lively 
debate lately 0Jl^l ■ One way to include such dephasing 
processes in the theory is provided by the Biittikker phe- 
nomenological model This approach is remarkably 
successful and its use became customary in the treat- 
ment of conductance fluctuation in chaotic dots. [T^] 
It introduces a ficticious voltage probe lead </>, through 
which there is no net current flow, but allows for elec- 
trons to randomize their phases at the reservoir <fi. As a 
result the dimensionless conductance reads 



9l p(E,X) = T 21 



T2<fE<i>\ 

T\d> + T; 
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where the arguments E and X are implicit to T. 

Temperature also affects the conductance in another 
(rather trivial) manner: the electrons flowing through a 
quantum dot are thermally distributed, yielding 



dE 9lP (E,X) 



dE 



(3) 



where / M = {1 + cxp[(_E — ^)/T]} -1 is the Fermi distri- 
bution function and /i is the chemical potential of the 
dot. It should be emphasized that in our notation g^ ac- 
counts solely for dephasing, while G is affected both by 
dephasing and the smearing of the Fermi surface. 

Let us consider the simplest statistical measures of 
P(G), namely the mean conductance (G) and its vari- 
ance var(G). In experiments (G) is obtained by varying 
fi and/or X, whereas in theory one takes a suitable en- 
semble averaging over g^. Actually, from Eq. (||) it is 
evident that (G) — (e 2 /h)(g^,). The inspection of the 
conductance autocorrelation function 



C{x) = (G(fi,X+)G((x,X-)) - {G(^X)f 



(4) 



where X^ = X ± x/2, reveals that the relation between 
var(G)[= C(0)] and var(g^) is obtained from M 



(kj C<j) {u,x) — 



2Tsinh 



(-) 

\2TJ 



(5) 



where c^iu, x) is the dimensionless conductance autocor- 
relation function defined by 



Cjfax) = (g^E+,X+) g4> (E-,X-)) - (g^ 



(G) 



with E ± = E±u/2. 

So far all theoretical studies [f^]8[|l6|-{il| aiming to de- 
scribe P{G) used the information theoretical approach. 
More specifically, one finds the ensemble of S matrices 
that maximizes the information entropy and fulfills the 
symmetries and other constraints of the physical system 
under consideration. This procedure is very amenable 
for the analytical calculation of -P(g^) (at fixed E and 
X) but has the limitation of lacking parametric corre- 
lations (neither for E nor for X) between members of 
the ensemble. This ingredient is of central importance 
in obtaining the variance and higher moments of G as 
indicated by Eq. (||). To circumvent this problem Ref. 
jl| introduced an heuristic procedure of smearing P(G), 
which underestimates var(G)^ = 2/var(G) i a = i. 

We use the Hamiltonian approach to the statistical S- 
matrix instead pc|| . Both frameworks are equivalent for 
the calculation of var(g^) |2l[], but not for var(G). The 
resonant S'-matrix is given by 



S(E,X) = t-2niW i 



1 



E - H(X) + iirWWi 



W (7) 



where H is taken as a member of the Gaussian or- 
thogonal (unitary) ensemble for systems where time- 
reversal symmetry is manifest (absent). For simplic- 
ity we take the case of N open channels at each lead. 
Since the H matrix, of dimension M , is statistically in- 
variant under orthogonal {(3 — 1) or unitary (/3 = 2) 
transformations, the statistical properties of S depend 
only on the mean resonance spacing A determined by H 
and on W , the coupling matrix elements between res- 
onances and channels. Those enter the theory through 
y c = tt 2 (W^W) cc /(MA) contained in the so called stick- 
ing coefficients P c = Ay c / (l + y c ) 2 . The later quantify the 
transmission through a given channels c, being maximal 
for P c = 1. By assuming the channels to be equivalent 
we can drop the index c. For open quantum dots the 
transmission is large and consequently it is assumed that 
P«l. In addition, we consider N^, open channels at the 
voltage probe lead, each of them with a sticking coeffi- 
cient p. The loss of phase coherence is modeled by the 
single parameter P^ = pN^, with S> 1 0]. The later 
can be converted in a dephasing width = AP^/27r, 
from which the dephasing time = h/T^ is extracted. 

In this approach the parametric correlations are auto- 
matically taken into account, but due to technical reasons 
it is very difficult to proceed analytically, unless N 3> 1 
and = 2 p9| . On the other hand, numerical simulations 
can be implemented straightforwardly. For each realiza- 
tion of H we invert the propagator and compute S for 
values close to the center of the band, E = 0, where the 
A is approximately constant. The dimension of H was 
fixed to be M — 200, taken as a compromise between 
having a reasonable wide energy window to work with 
and not slowing too much the computation. For each 
case under consideration we obtained good statistics for 
P(g<p) and c^(u!,x) with 10 4 realizations. 

We find that for the case of experimental interest, 
N = 1 and Pfy ^ 0, the numerically computed dimen- 
sionless autocorrelation function c$(u>, x) is quite similar 
to the one obtained in the semiclassical regime (N 3> 1) 



and P = lit 



namely 



var 



[1 + {x/X c ff + {uj/Tf 



(8) 



Our results, shown for the (3=1 case in Fig. [I] ((3 = 2 
gives essentially the same agreement), scale according to 



A 



r = r + -f = — 2NP + -f 



(3 2ir 



r 



(3 



(9) 



and X c = Kyj2NP + P^, where k is system specific and 
depends on the kind of parametric variation. It is worth 
mentioning that for P^ = there is additional work |24| ] 
showing that Eq. (j^) is a good approximation for any N. 
For simplicity we take P = 1 for the moment. 
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FIG. 1. Ratio between the dimensionless conductance au- 
tocorrelation function Cf,(x,uj) and its variance 0^(0,0) as a 
function of (a) x/X c with u) — and (b) uj/Y with x — 0. The 
solid line stands for the result of Eq. where as the numer- 
ical simulations for P$ — 1.0, 2.0, 3.0, and 4.0 are represented 
by the symbols x, □, o, and A respectively. 



The approximation for c$(u>, x) given by Eq. 
for an analytical evaluation of C[x) in Eq. (||) 



C(x) 



eT 



2var(g^,) 



x z + 



1 hi & 2 + 1 + n2^T/Tf 
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FIG. 2. Ratio between var(G) and var(g^) (in units of 
e 4 /h 2 ) as a function of T/r. The sum in Eq. (|lo|) is of slow 
monotonic convergence. The thick solid line gives the con- 
verged result, while the two thin ones stand for evaluation of 
the first two and five terms of the sum. The dotted line is the 
approximation given by Eq. (pd[). 



where x — x/X c . For C(0) = var(G), the above expres- 
sion is nicely approximated (within 15%) by 



2\ 2 

e \ var 



Var(G)= lXy 1 + 2T/F 



(11) 



as shown in Fig. In Eq. ( pd| ) the dependence of 
var(G) on /3 is implicit in both var(<7^) and T^. Figure 
U shows that these considerations reconcile the theory 
with the experimental data. The information theoretical 
approach underestimates the ratio var(G)^ = 2/var(G) ) g = i 
because it lacks the temperature correction given by 
(1 + 2T/r (3=2 )/(l + 2T/r i a=i). Notice that we do not 
introduce any additional fitting parameter in our theory. 




FIG. 3. Ratio var(G)/3=2/var(G),3 = i as a function of P$. 
The dashed line stands for the information theoretical ap- 
proach interpolation formula, while the circles (•) are the ex- 
perimental data, and the solid line is our result. 

We now switch our attention to the question of how to 
determine the dephasing rate in open chaotic quantum 
dots. There are three main proposed strategies p5| , pT| . 
Let us start addressing the one based on the weak lo- 
calization peak. As shown (G) = (e 2 /h)(g^), allows one 
to read the average dimensionless conductance directly 
from the empirical data. In turn, provided that the leads 
are ideal (P — 1), the weak localization peak, defined as 



&9 = (.90>/3=2 - (94>)/3=i 



(12) 



is in direct relation to T^. 11 The problem is that in ac- 
tual experiments P < 1. Thus, 5g is a function not only 
of but of P as well. An inspection of Fig. |], obtained 
from our simulations for N = 1, shows that by fixing Sg 
(as obtained from the experiment) and decreasing P by 
a small factor always increases Pa. The effect becomes 
small for P^ 1, but is rather large for P^ < 1. In this 
situation, reducing P by 5% decreases P^ by as much 
as 100%. The dephasing time increases accordingly. 



3 



Hence, Sg does not uniquely fixes P^. This ambiguity can 
by eliminated by using the experimental (g$) for j3 = 1 
and 2 to fix both coefficients P and P^. The data from 
Ref. pi indicate that the correction to P^ is significant. 



of T0 from the data, and propose an alternative solution. 
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FIG. 4. Contour plot of the weak localization peak 8g as a 
function of the coefficients P,/, and P for N = 1. 

There are two other methods to extract from ex- 
periments dealing with chaotic quantum dots. Both are 
related and rely on the study of the parametric depen- 
dence of the conductance. Based on a semiclassical argu- 
ment [ p5| , it was proposed that the study of the conduc- 
tance autocorrelation as a function of an external mag- 
netic field, C(x = B/B c ), has a simple dependence on 
r , namely (B c /$ ) 2 = n(2NP + (Our numeri- 

cal results support this relation, as depicted in Fig. [I].) 
By measuring B C (T) one can thus obtain t^,(T). The 
problem here is that C(x) changes its functional depen- 
dence when going from the T ^ A to the T <C A limit 
p9[ , which can jeopardize the determination of B c for 
T « A. Alternatively, the width of the Lorentzian dip 
of the average conductance around B = 0, can also be 
used ( 94> (B)) = (^>/3=2 - 6g/[l + (B/B c ) 2 }. Both meth- 
ods were recently shown to give consistent results with 
the weak localization one, at least for T > A fLlf . This 
is in apparent contradiction with our claim that there 
is an ambiguity in the weak localization peak method. 
However, differences are only expected for small values 
of P^, where the parametric methods were not employed. 
Moreover, since N = 1 and P w 1 give B c oc 2 + P^ , 
the latter become evidently inaccurate for P^ <C 1. 

In conclusion, we presented a detailed statistical study 
of conductance fluctuations in chaotic quantum dots. We 
solved the only serious discrepancy between theory and 
experiment, giving a stronger support to the statisti- 
cal approach incorporating dephasing. In addition, we 
pointed out some problems in the quantitative assertion 
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